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Abstract For many spin systems with constant isotropic antiferromagnetic next- 
neighbour Heisenberg coupling the minimal energies E„iin{S) form a rotational 
band, i. e. depend approximately quadratically on the total spin quantum number 
S, a property which is also known as Lande interval rule. However, we find that 
for certain coupling topologies, including recently synthesised icosidodecahedral 
structures this rule is violated for high total spins. Instead the minimal energies are 
a linear function of total spin. This anomaly results in a corresponding jump of the 
magnetisation curve which otherwise would be a regular staircase. 
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1 Introduction 

It appears that for spin systems with constant isotropic antiferromagnetic next- 
neighbour Heisenberg exchange the minimal energy Emin (5*) for given total spin 
quantum number S is typically a strictly convex function of S. For many spin 
topologies like rings, cubes, icosahedra etc. this function is very close to a parabola 
[[l|]. For certain systems this behaviour has been explained with the help of the un- 
derlying sublattice structure [^. Experimentally this property has been described 
as "following the Lande interval rule" In the classical limit, where the 

single-spin quantum number s goes to infinity, the function Emin{S) is even an 
exact parabola if the system possesses co-planar ground states 

However, we find that for certain coupling topologies, including the cubocta- 
hedron and the icosidodecadron [Kfl, this rule is violated for high total spins. More 
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precisely, for the icosidodecadron the last four points of the graph of Emin versus 
S, i. e. the points with S = Smax to 5 = Smax — 3, lie on a straight line 

EminiS) = eOJs^ _ 6Js(30s - S) . (1) 

An analogous statement holds for the last three points of the corresponding graph 
for the cuboctahedron. These findings are based on numerical calculations of the 
minimal energies for several s both for the icosidodecahedron as well as for the 
cuboctahedron. For both systems, additionally, we have a rigorous proof of the 
high spin anomaly for the case of s — 1/2. This proof rests on an inequality 
which says that all points of the graph of E,nin versus S lie above or on the line 
connecting the last two points ("bounding line"). The proof can be easily applied 
to a wide class of spin systems, e.g. to two-dimensional spin arrays. 

The observed anomaly - linear instead of parabolic dependence - results in 
a corresponding jump of the magnetisation curve Ai versus B. In contrast, for 
systems which obey the Lande interval rule the magnetisation curve at very low 
temperatures is a staircase with equal steps up to the highest magnetisation. 

The anomaly could indeed be observed in magnetisation measurements of the 
so-called Keplerate structure {Mo72Fe3o} which is a recently synthezised mag- 
netic molecule where 30 Fe'^+ paramagnetic ions (spins s = 5 /2) occupy the sites 
of a perfect icosidodecahedron and interact via isotropic, next-neighbour antifer- 



romagnetic exchange [|9p. Unfortunately, the magnetisation measurements 1 1 1 
performed so far suffer from too high temperatures which smear out the anomaly. 

Nevertheless, it may be possible to observe truely giant magnetisation jumps in 
certain two-dimensional spin systems which possess a suitable coupling topology. 
In such systems the magnetisation jump can be of the same order as the number 
of spins, i.e. the jump remains finite - or is macroscopic - in the thermodynamic 
limit N ^ oo. 

The article is organized as follows. In section || we introduce basic definitions 
and explain how the results have been obtained. In section^ the high spin anomaly 
is discussed and proven for the case of s = 1/2. We provide an outlook in section 

i 



2 Definitions and Numerics 

The Heisenberg Hamilton operator of the investigated spin systems is 

^=IY1 -^"^ • + 9f^BBS, (2) 

(u,v) 

= ^ J2 • s(w) +.g^BSS;^ , S^ = '^s^{u), 

iu,v)er « 

where the exchange parameters J„ „ are considered as the components of a sym- 
metric matrix J, i.e. every bond is taken into account twice. In particular, we as- 
sume J„ „ S {J, 0} and J > which corresponds to antiferromagnetic coupling. 
In equation (0), g is the spectroscopic splitting factor and jj, b the Bohr magneton. 
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The vector operators s (u) are the spin operators (in units of K) of the individual 
N paramagnetic ions with constant spin quantum number s. Because the matrix 
J couples only next-neighbours (see Fig. |l]) the second sum in (^ runs over the 
set r of all next-neighbour pairs {u, v) of spins of a single molecule at sites u 
and V. r can be regarded as the set of "edges" of the corresponding undirected 
graph describing the coupling scheme of the molecule. The "vertices" of the graph 
correspond to the spin sites 1, 2, . . . , iV. For each spin site u let r{u) denote the 
set of neighbours of u. Throughout this article we will assume that the number of 
neighbours per site is constant, say = j- The "distance" between two spin 

sites u and v will be the minimal number of edges connecting u and v (similar to 
the Manhattan distance). 




As mentioned already in the introduction the anomaly was found numerically. 
For this purpose the Hamilton matrix had to be diagonalized. The total matrix is a 
huge object of dimension (2s + l)^x (2s + l)^ which must be block-diagonalized 
in advance. Using that the Hamilton operator commutes with Sz, the Ising product 
states which are a natural basis can be grouped according to the quantum numbers 
M, thereby dividing the Hilbert space into orthogonal subspaces 7i{M). A further 
reduction of dimension is achieved if the symmetries of the spin array are ex- 
ploited. The icosidodecahedron for instance shows a tenfold shift symmetry lead- 
ing to Hilbert subspaces 7Y(il/, k) with k — 0, . . . , 9. Within these subspaces a 
Lanczos procedure was applied in order to obtain the respective minimal energies. 



3 High spin anomaly 

3.1 Observations 

The resulting minimal energies Emin{S) are shown by dashes on the l.h.s. of Fig. ^ 
for the isosidodecahedron and on the l.h.s. of Fig. ^ for the cuboctahedron. The 
straight lines denote the bounding lines, which connect the highest four levels in 
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the case of the isosidodecahedron and the highest three in the case of the cubocta- 
hedron. At T = this behavior leads to jumps of the magnetisation A4 

X = -| tr [gfiBS.e-^'" } , Z ^ U [e-^'i } . (3) 

Due to the effect that the states he exactly on the bounding line in the graph of 
Emin versus S they "take over" for the new total ground state at the same value of 
the magnetic field, therefore the magnetisation immediately jumps to the highest 
value. The jumps are marked by arrows in the magnetisation curves of the isosido- 
decahedron (r.h.s. of Fig. ^ and the cubeoctahedron (rh.s. of Fig. ||). 





2 4 6 8 10 12 U 
S 



5 10 1.5 2 2.5 3 3 5 



Fig. 2 Icosidodecahedron: L.h.s. - minimal energy levels Emin{S) as a function of total 
spin S. R.h.s. - magnetisation curve at T = 0. 





5 10 15 2 2 5 3 3 5 



Fig. 3 Cuboctahedron: L.h.s. - minimal energy levels Emin{S) as a function of total spin 
S. R.h.s. - magnetisation curve at T = 0. 



For systems which follow the Lande interval rule, i.e. where Emm{S) is a 
parabolic function of S, the corresponding magnetisation curve would consist of 
equal steps. 



3.2 Idea of the proof 

A necessary condition for the anomaly is certainly that the minimal energy in the 
one-magnon space is degenerate. Therefore, localized one-magnon states can be 
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constructed which are also of minimal energy. When placing a second localized 
one-magnon eigenstate on the spin array there will be a chance that it does not 
interact with the first one if a large enough separation can be achieved. This new 
two-magnon state is likely the state of minimal energy in the two-magnon Hilbert 
space because for antiferromagnetic interaction two-magnon bound states do not 
exist (at least for s = 1/2). This procedure can be continued until no further 
independent magnon can be placed on the spin array. In a sense the system behaves 
as if it consists of non-interacting bosons which, up to a limiting number, can 
condense into a single-particle ground state. 

In more mathematical terms: In order to prove the high-spin anomaly we first 
show an inequality which says that all points {S, Emin{S)) lie above or on the 
Une connecting the last two points. This inequality holds for s = 1/2 and all sys- 
tems with constant antiferromagnetic exchange parameter and a constant number 
of neighbours for each spin site. For specific systems as those mentioned above 
what remains to be done is to construct particular states which exactly assume the 
values of Emin corresponding to the points lying on the bounding Une, then these 
states are automatically states of minimal energy. 

Note that the high spin anomaly does not contradict the strict convexity of 
the graph of Emin versus S in the classical limit, since in the limit .s oo the 
interval where the anomaly occurs, e.g. S = Smax — 3 to 5 = Smax, becomes an 
infinitesimally small fraction of the total spin range. 

We set J = 1 throughout this section. 



3.3 Bounding line for s = 1/2 

Let Jia denote the eigenspace of Sz with eigenvalue M ~ N/2 — a, a = 

0, 1, . . . , iV. It has the dimension dim(7ia) = ( ^ 1 • orthonormal basis of 

Ha is given by the product states denoted by \n\,. . . , Ua) with 1 < m < n2 < 
■ ■ ■ < ria < N where the rij denote the sites with flipped spin m = — 1/2. A state 
of this form will be called isolated iff {rii, rij) ^ F for all 1 < « < j < a. In 
other words, the flipped sites of an isolated state must not be neighbours according 
to the coupling scheme. will denote the subspace of Ha spanned by isolated 
states. 

We will embed the Hilbert space of the spin system into some sort of Fock 
space for magnons. More precisely, let Ba{Hi) be the totally symmetric (i. e. bosonic) 
subspace of 0^=1 ^Hi.JfAi : Hi — > Hi is a linear operator, Ba{Ai) will 

denote the restriction of^i^J^^i-'-^iJ^H \- l<g) ■ ■ ■ (g) 1(g) Ai onto Ba{Hi). 

An orthonormal basis of Ba{Hi) is given by the bosonic states 

S\ni)^---^ \na), (4) 

where 1 < rii < n2 < ■ ■ ■ < ria < N and the "synmietrisator" S denotes the sum 
over all a! permutations of the product state divided by \/aI. The Unear extension 
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of the map \ ni, . . . ,na) i-^ 5 | ni ) ® • • • \na) defines an isometric embedding 

Ja-.Ha^ BaiHl) . (5) 

Let Ha denote the restriction of the Hamilton operator (with zero magnetic 
field) onto Tia and 

Ha = JlBa{Hl)Ja- (6) 

We will show the following 
Proposition 1 Ha = ^^Nj + Ha if restricted to the subspace Tia°- 

j is the number of neighbours, which is assumed to be constant for each spin site. 

Proof: Let | ni, na ) be an arbitrary isolated basis state and split s{u) ■ 
s{v) into 

s{u) ■ s{v) = s^{u)s^{v) + ^{s + {u)s-{v) + s-{u)s + {v)) , (7) 
analogously Ha = H' + H ' and Ha = Ha + Ha ■ First, let us consider 

H';\m,...,na)^^ E is + {u)s-{v) + s-{u)s + {v))\n,,...,na)(S) 

= M E Sort I mi,..., na ) H h ^ Sort | rii, . . . , Wa ) 

Here Sort denotes the procedure which re-arranges a list of integers into its non- 
decreasing order. Note that further summation constraints of the form rni ^ n2, 
. . . ,na etc. would be superfluous since \ni, . . . ,na) was assumed to be isolated. 
Now consider 

Ha" \ ni, . . . ,na) ^ JlBaiH'DSlm) (g) ■ ■ ■ (g) \na) (9) 
I mi ) (g) • • • (g) \ na) 



\ mier(ni) 



+ • • • + E 5 I m ) (g • • • (g I ma ) 

'"aer(n„) / 

\\ E Sort I mi,..., na ) H h E Sort | m 



, ■ ■ ■ , I'^a I 



4 

= Hl\n^,...,na) . (10) 

Now we turn to H'^. Recall that there is a total number of L = ^ links between 
different sites. For a given basis state | ni, . . . , na ) we write 

L = L++ + L+_+L__, (11) 
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where denotes the number of links between two rn = +l/2-sites, etc. Hence 
for isolated states = 0. Each basis state is an eigenstate of H'^ with eigen- 
value 

i(L++-L+_+i__) = i(L-2L+_). (12) 

For isolated states i_| = ja, but in general i_| = ja — 2i__ since each 

link "deletes" two H — links of a corresponding isolated state. Hence 

H'a\m,...,na) = j - 2ja + 4L__^ \ ni,...,na) (13) 

= J (^-2ja^ |ni,...,na). (14) 

Similarly one can show that 

Ha \ni,...,na) = -{— - 2j j | ni, . . . , ti^ )■ (15) 

From ( 15 ), (p4|), and (10) the proposition follows immediately. ■ 
If we drop the condition that | ni, . . . , ) is isolated we ought to slightly 
modify our calculations. First, we would have to introduce extra summation con- 
straints of the form mi 7^ ri2 , . . . , Tia etc. in (^ in order to be sure that the resulting 
states lie in Ha- Although Ba{Hi) in general will produce some unphysical states 
with two magnons locahzed at the same site, these states will be annihilated by 
Jl. Hence again H'^lm, . . . ,na) ^ Ha | m, . . . , ). 

For H[^ and Ha the situation is different. The foregoing arguments show that 
the difference between the left and the right hand side of proposition |l] is some op- 
erator with eigenstates | ni, . . . , ) and corresponding eigenvalues L which 

are > 0. This shows that proposition |l]generaUzes to 

Proposition 2 Ha> ^Nj + Ha- 

Now let Ea denote the smallest energy eigenvalue of Ha- Note that (5* = 
N/2 — a) > Ea, since the energy eigenvalues for given total spin quantum num- 
ber S are assumed within each subspace of magnetic quantum number M — 
—S,-.-,S- We expect that Emin{S — N/2 — a) — Ea holds generally for the 
spin systems under consideration, this has been proven only for so-called bi-nartite 
systems [ 12, 13|], but numerically shown to hold for much more systems [[l4{|. 



Analogously we define Ea for Ha- Since for bosons the ground state energy 
is additive, aEi will be the smallest energy eigenvalue of B{Hi). We further con- 
clude > aEi since Ea - {^\Ha\^) = {Ja<^\B{Hi)\Ja'P) if Ha\<P) = 

i?o I ^ ) ■ Together with proposition |l] this implies 

Proposition 3 Ea>^Nj + aEi . 

This inequality says that the minimal energies Ea, resp. Emin{S — N/2 — a), 
he above or on the "bounding line" £{a) = ^-g^Nj + aEi. 
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3.4 Ground states of independent magnons 

According to what has been said above in order to rigorously prove the high spin 
anomaly it suffices to construct states which assume the energy values of the 
bounding line £{a) for certain values of a > 1. By the results of the previous 
subsection it is clear that these energy values must be minimal and the states must 
be eigenstates of Ha in the case s = 1/2. Actually we conjecture that these states 
are also minimal energy states of Ha for arbitrary spin, which conjecture is numer- 
ically supported for all cases where we have calculated Ea, but we cannot prove it 
at the moment. Nevertheless, we will assume an arbitrary spin s in this subsection. 

We first consider the case of the icosidodecahedron. Let a — 1. Recall that 
a general state in Ha is of the form X]^=i where n denotes the spin site 

where the magnetic quantum number is decreased by 1 . The eigenvalues of Hi are 
of the form 

Eo. = ^Njs^ + {j^-j)s, (16) 

where ja,a — I, . . . , N aie the eigenvalues of the coupling matrix J. In our case, 
N = 30, j — 4, and the minimal eigenvalue ja is —2, hence 

El = 60s^ - 6s. (17) 

The corresponding eigenspace of Ha is ten-fold degenerate. It is possible to find 
linear superpositions which are states of minimal energy and have some intu- 
itive geometric interpretation as localized one-magnon states corresponding to 
even subrings of the icosidodecahedron. These states have alternating amplitudes 
c„ = ±1 for sites n of the subring and vanishing amplitudes for the remaining 
sites. The smallest even subrings generating such states are the "8-loops" circum- 
scribing two adjacent pentagons, e. g. (1,2,3,4,9,17,18,10) according to the 
numbering of sites in Fig. [l| Other even subrings are the "equators" with 10 sites 
or the "curly equators" with 12 sites which need not be further considered here. 

Now let a = 2. If a two-magnon ground state lies on the bounding line £{a), 
as it is suggested by numerical diagonalization, we would have 

E2 = 60s^ - 12s. (18) 

In fact, this energy is assumed by the following state: Consider two 8-loops Li , ^2 
with a distance of 2, e. g. Li = (1,2,3,4,9,17,18,10) and 
L2 = (12, 13, 22, 28, 29, 30, 26, 21) according to Fig.[l|. e^^ , € ii and S^^ , 712 G 
L2 denote the amplitudes which define the one-magnon ground-states described 
above. Then a two-magnon ground-state with the energy can be defined by 

<^2 = eni'5„Jni,n2). (19) 

niGii,n2Si2 

Since this state lies entirely in Hf" it can be considered as a ground-state of two 
non-interacting magnons. 
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Table 1 Definition of an 3^h-symmetric three-magnon ground-state by assignment of am- 



plitudes to representative triple states. 

1 ni , n2 , 713 ) Length of orbit Amplitude 


1 1,3, 14) 


60 


1 


|1,3,15> 


120 


-1 


1 1,3,22) 


120 


-1 


1,3,23) 


120 


1 


1,3,28) 


120 


1 


1,3,29) 


60 


-2 


1,7,15) 


120 


1 


1,7,18) 


60 


1 


1,7, 23) 


120 


-1 


1,7, 24) 


120 


-1 


1,7, 29) 


60 


2 


1 1,8, 21) 


120 


-1 


1 1,8,25) 


30 


2 


1 1,8, 26) 


120 


-1 


1 1,8,27) 


120 


2 


1,8,28) 


30 


-2 


1,13, 16) 


20 


-2 


1,13,24) 


60 


-1 


1,13,30) 


60 


2 


1,14, 30) 


20 


-1 



Unfortunately, an analogous construction of three mutually isolated one-magnon 
states is no longer possible for a = 3. Here we have to determine an appropri- 
ate state by numerical diagonalization. One possible state of three independent 
magnons is a state which is completely symmetric under the action of the sym- 
metry group of the icosidodecahedron, i.e. the icosahedral group with reflections 
yh of order 120. Hence it will suffice to define this state by assigning an ampli- 
tude to only one triple of sites within each orbit of the symmetry group. The other 
triples obtained by applying symmetry operations g ^ yn lo each site will have, 
by definition, the same aplitude. The complete definition of this state (without nor- 
malization) can be found in Table |l]. The calculation of the corresponding energy 
E2 = 6Qs^ — 18s can be done by a computer algebra software. Also this state Ues 
entirely within TC^^". Thus we have obtained a rigorous proof of the anomaly also 
for the case a — 3 and s = 1/2. 

The case of the cuboctahedron is largely analogous, up to the fact that here 
we have only one point of anomaly for a = 2. The corresponding two-magnon 
ground state can be constructed by using two separated 4-loops, e. g. (1, 2, 3, 4) 
and (9,10,ll,12)inFig.|l](r.h.s.). 
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3.5 Generalization to the XXZ-model 

The above proof holds also for the more general Hamiltonian of the XXZ-model 

H=^^ J2 {^lz{u)s,{v) + s^{u)s^{v) + Sy{u)sy{v)^ , (20) 

iu,v)er 

with A > Q. Since the total spin S is no longer a good quantum number, the 
minimal energies E„iin have to be considered as a function of the total magnetic 
quantum number M instead. For the existence of the bounding line and the cor- 
responding magnetisation jump this aspect is irrelevant. The only change in the 
proof is a multiplication of H'^ and Ha by A, which does not change the argu- 



mentation. Also the construction of eigenstates, as carried out in subsection 3.4, is 
not altered by the anisotropy A in (EG), since these states are isolated. 



4 Outlook 

The shown proof offers a method to create spin arrays which by construction sup- 
port a finite number of independent magnons. The basic idea is to design a unit cell 
which can host a localized one-magnon state, that is an eigenstate of the Hamilto- 
nian. Triangles play a key role in the construction of such cells because they help 
to prevent localized magnons from escaping. The total spin array is then obtained 
by properly linking serveral unit cells. Fig. ^ shows an example. The unit cell is 
one quarter of the structure. It can host a single magnon 

llmagnon) = i(|l)- |2)+ |3)- |4)) , (21) 

which is an eigenstate of the Hamiltonian with minimal energy in the one-magnon 
space. One easily notices that in total four localized independent magnons fit into 
the structure. In general it might be possible that more independent magnons, like 
in the case of the icosidodecahedron, can occupy the spin array. For the example 
of Fig. this is not the case. 

The latter example offers the perspective of observing truely giant magnetisa- 
tion jumps in two-dimensional spin systems. The number of independent magnons 
which can be placed on the lattice is proportional to the number of spins itself - 
N/6 is the example of Fig. ^- and thus a macroscopic quantity. This will be the 
subject of a forthcoming pubhcation. 

Acknowledgements J. Richter and J. Schulenburg thank the Deutsche Forschungsgemein- 
schaft for support (project Ri 615/10-1). 
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